arXiv: 1507.05077vl [cond-mat.mes-hall] 17 Jul 2015 


Magnetotransport in Dirac metals: chiral magnetic effect and quantum oscillations 

Gustavo M. Monteiro,^ Alexander G. Abanov/’^ and Dmitri E. Kharzeev^’^’^ 

^Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794, USA 
^Simons Center for Geometry and Physics, Stony Brook University, Stony Brook, NY 11794, USA 
^Department of Physics, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
^RIKEN-BNL Research Center, Upton, New York 11973-5000, USA 

Dirac metals are characterized by the linear dispersion of fermionic quasi-particles, with the Dirac 
point hidden inside a Fermi surface. We study the magnetotransport in these materials using chiral 
kinetic theory to describe within the same framework both the negative magnetoresistance caused 
by chiral magnetic effect and quantum oscillations in the magnetoresistance due to the existence of 
the Fermi surface. We discuss the relevance of obtained results to recent measurements on Cd3As2. 


The discovery of Dirac semimetals w has enabled 
the experimental studies of 3-dimensional materials with 
chiral quasiparticles. In comparison to 2-dimensional 
graphene, the access to three spatial dimensions allows 
one to study phenomena such as chiral anomaly [MS] and 
chiral magnetic effect (GME) [9]. In particular, the lat¬ 
ter refers to the generation of electric current induced by 
chirality imbalance in the presence of the magnetic held; 
see unmn for reviews. 

Dirac semimetals together with Weyl semimetals are 
representatives of three-dimensional chiral materials. 
They are both characterized by the existence of band¬ 
touching points. In Weyl semimetals, the time reversal 
symmetry is spontaneously broken and each Dirac node 
splits into two disjoint Weyl points. In a Dirac semimetal, 
each Dirac node has zero Chern number, reflecting the co¬ 
existence of two Weyl points of opposite chiralities at the 
same point in the Brillouin zone. This corresponds to an 
emergent Z 2 -symmetry relating the two states with dif¬ 
ferent chiralities. In parallel electric and magnetic helds, 
the degeneracy between the states with opposite chirali¬ 
ties gets broken due to the chiral anomaly and the differ¬ 
ence between the Fermi energies of left- and right-handed 
fermions can be described by the chiral chemical poten¬ 
tial This difference generates a nonvanishing chiral 
magnetic current of the form Jcme = PbB/{2 'k'^). 

Because of the chiral anomaly, /j,^ ^ E ■ B. Conse¬ 
quently, the CME conductivity acquires a positive term 
proportional to and the magnetoresistance (MR) 
becomes negative m- Such a behavior signaling the 
presence of CME has been observed recently in Dirac 
semimetals ZrTes [13] and NaaBi [14]. Negative MR has 
also been observed [Hide] in TaAs, a candidate for a 
Weyl semimetal. On the other hand, the previous stud¬ 


ies of magnetotransport in another candidate for a Dirac 
semimetal Cd 3 As 2 [niiH] revealed a more complicated 
pattern, with strong oscillations of MR. 

Quantum oscillations in MR signal the presence of 
large Fermi surface. The material represents a Dirac 
metal rather than a semimetal with the Dirac point hid¬ 
den inside a Fermi surface. In this Letter we develop 
the theory of magnetotransport in Dirac metals that 
describes an interplay between the CME and quantum 
Shubnikov-de Haas (SdH) oscillations. Our analysis can 
be trivially extended to Weyl metals if one assumes only 
one chirality per Dirac node. 

Although Dirac metals are characterized by linear dis¬ 
persion of quasiparticles e{k) = hvF\k\*, the assumption 
of linear spectrum is absent in this Letter. However, 
we assume that the quasiparticles can be described by 
the Fermi liquid theory. Also, we restrict ourselves to 
isotropic system. The chemical potential or Fermi en¬ 
ergy define the size of Fermi surface ep = sikp), where 
the Fermi momentum is related to the density of conduc¬ 
tion electrons (per Dirac point and chirality) by standard 
formula kp = (37r^ne)^/^. 

Our treatment will be based on the semiclassical ap¬ 
proximation, valid when magnetic field is weak enough 
so that a large number of Landau levels is filled. Intro¬ 
ducing the magnetic length ip = this condition 

amounts to ^kpi^ 1. In this limit, we can associate 
a trajectory to each electron quasiparticle in the vicinity 
of the Fermi surface.t We expect to see pronounced SdH 


* We assumed the Dirac point to be at fc = 0. 

^ In this limit one can still think about Fermi sphere albeit strati¬ 
fied into Landau level “cylinders”, see Figure [T| 
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oscillations when the temperature is much smaller than 
the energy gap between Landau levels or, equivalently, 
when the thermal de Broglie wavelength is much bigger 
than the Larmor radius 


A 




FIG. 1: Semiclassical picture of Fermi surface for right and 
left chiral modes is shown in the presence of the magnetic 
field. For the linear dispersion, the support of Hamilto¬ 
nian eigenstates is a collection of cylinders corresponding to 
eigenvalues en(fcz) = hvpy/kl ~\- k\, where k\ = 2eBn with 
n = 0,1, 2,.... The state with n = 0 is chiral and exists only 
for Pz > 0 (parallel to B) for the right chirality and for Pz < 0 
for the left one. 

At low temperatures, the impurity scattering is the 
leading contribution to conductivity tensor. This intro¬ 
duces another scale into the problem - the scattering 
rate. The system is called clean when the quasiparticle 
performs many cyclotron orbits before colliding or, equiv¬ 
alently, when the mean-free-path is much larger than the 
cyclotron radius, vpT 3> We restrict ourselves to 

single-impurity scattering approximation and neglect in¬ 
terference and localization. This approximation is valid 
when the density of impurities is low. 

Thus the regime of interest in this work is defined by 

1 <C ‘C kpvpT , (1) 

where for the linear spectrum the last term correspond 
to epT/h. 

Dirac points in 3 dimensions correspond locally to 
monopole solutions of the Berry curvature in momen¬ 


tum space m, which give rise to chiral anomaly effects 
in kinetic theory [20l [21] ; see [22] for application to Weyl 
semimetals. The Berry curvature is obtained from the 
Bloch functions for the valence band electron quasiparti¬ 
cles \uk) and for an isotropic system is given by ^5] : 


k 

r2(fc) = iVfc X (wfclVfeMfc) = . (2) 

The Berry curvature ([^ has an opposite sign for chiral¬ 
ities X = if and k denotes a unit vector in the direction 
of k. We refer to the modifications of kinetic theory by 
Berry curvature as to the chiral kinetic theory iiniEi]. 
In particular, the Berry curvature modifies the expression 
for the current density. Given a dispersion relation e(fc), 
we introduce the group velocity vector = ^Vfce(fc), 
so that the current density can be expressed as [24U27] : 


j = -2e / / 


Vk p i^k ' ^ ^ ^ 

n n 


d^k 
(27r)3 ’ 


( 3 ) 

where /(a;, k,t) is the distribution function and the inte¬ 
gral is performed over the first Brillouin zone (BZ). The 
overall factor of 2 accounts for spin projections. 

The distribution function is obtained by solving the 
Boltzmann equation. Since we are interested in linear 
response, we expand f{x,k,t) around the equilibrium 
(Fermi-Dirac) distribution function fo{e): 


f{x,k,t) = /o(e) + e 


dfo 

de 


E-g + 0{E^). 


( 4 ) 


Having g{x,k,t) from the linearized Boltzmann equa¬ 
tion (see below) and substituting the ansatz Q into (§, 
the conductivity tensor reads: 



Here and in the following all the expressions refer to a 
single Dirac node. The assumption that the Dirac points 
can be treated independently is valid when they are far 
apart in the Brillouin zone^, so that the quasiparticle 


I We measure temperature in energy units. 


S In comparison to the Fermi momentum of each disjoint piece of 
Fermi surface. 
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scattering from one Dirac cone to the other requires a 
large momentum transfer. The last term in Eq. ([^ van¬ 
ishes for isotropic dispersion relations. The last equality 
is obtained with the use of assuming that the system 
is isotropic and that the integral is dominated by a vicin¬ 
ity to the Fermi surface due to the factor dfo/ds. We also 
introduced the small parameter (cf. §) a = i/i2kH b) 
and considered that the magnetic field is along the z- 
direction. 

For Dirac metals, the Z 2 -symmetry holds at low en¬ 
ergies and interaction terms that break this symmetry 
are sub-leading in comparison to the chirality-preserving 
ones. In this limit, the Boltzmann equations for differ¬ 
ent chiralities decouple and the collision integral accounts 
only for intra-chirality scattering. Given the transition 
rate Wk'^k from an initial state k' to a final state fe, the 
collision integral can be written as: 

m= [ (7) 

Jki, 

Here we assumed the elastic scattering probability to 
be invariant under time reversal, i.e. = Wk^k'. In 

addition to that, we have used [/' (1 — /) — / (1 — f')] = 
f{k') — f{k) and denoted the modification of the phase- 
space volume element due to a non-vanishing Berry cur¬ 
vature as 


T = l+^B-n{k) = l + xCkCose, (8) 

where 9 is the angle between k and magnetic field. 

Using the equations of motions from [MIET] , the Boltz¬ 
mann equation for g{t, k) in the linearized regime be¬ 
comes: 


Tidt+iuj)--{vkxB)-Vk 
e. f d^U 

= v, + -{v,^n)B + J 


a = 


(9) 


{T'wk'^kJ) [a' -a] ■ 


In Eq. (|^, we have assumed that the system is uni¬ 
form and the electric field oscillates with the frequency 
w, i.e., E = It is straightforward to observe 

that this equation does not admit a stationary solution 
when w = 0^. This is the manifestation of the chiral 
anomaly in kinetic theory - the constant parallel electric 
and magnetic field continue to pump chirality into the 


^ This can be seen by integrating (j^ over the solid angle. 


system. However, a stationary solution does exist in the 
presence of a chirality relaxation mechanism. 

To determine Wk'-^k^ we assume that the elastic scat¬ 
tering occurs on weak, dilute, and point-like impurities. 
We thus model the single-impurity scattering by 

Wk'~^k = „ . ^ {1 + k' ■ ie) S{e - e'), (10) 

2iy{£)T{e) 


where is the density of states - in the absence of 
magnetic field - at the energy e. We assumed that the 
scattering is elastic and averaged over impurity positions. 
All microscopic details are absorbed into the transport 
scattering time r. We remark here that although we 
focused on the small wave vector limit, the scattering 


rate from Eq. (10) is not isotropic. This is because the 


Weyl-particle spins are always polarized along their mo¬ 
menta, producing a universal factor (1 -|- • fc), which 

suppresses the backscattering of particles by impurities. 
For example, for massless Dirac quasiparticles one can 
find at leading order in the partial-wave expansion of 
scattering amplitude**: 


2vp 


1 _ 


sin^ <5i . 


The scattering phase (5i in the general case should also 
depend on the magnitude of magnetic field since the 
screening of the impurity potential might be modified 
by B. Since B = Bz, the azimuthal symmetry along 
the z-direction allows us to find solutions to that are 
independent of (j). Solving for gz{k,9), we find 


(k 9) = + 1-Cf Vk cos 9 

^ ’ iuj + g 1-1- xCfc cos 9 iuj -\-\/t ' 


( 11 ) 


where ry —> -1-0 in the absence of chirality flipping and 
will be replaced by 1 /r„ if the chirality flipping processes 
are taken into account. 

The phase space factor (|^ takes into account the redis¬ 
tribution of the density of states along the Fermi surface 
in weak magnetic fields, i.e., accumulation of states at 
the south/north pole for left/right chirality, respectively 


** Although the magnetic field breaks the 3D rotation invariance, 
the assumption of adiabatic evolution allows us to write the 
eigenbasis in terms of Bloch functions or plane waves. The ef¬ 
fect of magnetic field is absorbed into the trajectory in k-space 
and in the measure. A solution of the Dirac scattering prob¬ 
lem can be found, e.g., in m and gives for scattering amplitude 

A(k'-k) = i -1) [Pi(fc'-fc)-bD_i(fc'-fc)] .. 
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(see Fig. [^. However, it does not take into account the 
discreteness of Landau levels crucial for SdH oscillations. 
The discreteness of Landau level can be included through 
the Bohr-Sommerfeld quantization condition: 


6 \ 27r7Z 

1 + ■ B) z ■ k X dk = . (12) 

n / ^ B 


Eq. (12) conies from the non-trivial Poisson brackets 


between the coordinates of k and 7 denotes the curves of 
constant energy. Taking into account the discreteness of 
Landau levels into Eq. ([^, the conductivity per chirality 
becomes: 


o-(x) = _ 

^ ZZ 






n =0 


d(cos 6 ') /c^^Ufe(cos 6 » + xCk) 


-1 

c 2 


X i 5 f n — 


1 — cos"^ 9 — 2xCk cos I 

4Cfe 


gz{k,e). (13) 


Since the argument of the delta function has no real 
roots when n G Z_, we can consider the sum starting 
from n = —00 and use the Poisson summation formula. 
Thus, 


g .( x ) _ Q -( O ) 
^ZZ ^ zz 


2E 


(j(0 cos(— + - 


where we used dimensionless magnetic field 

Qf = (,k\k=kF = 


2kpl 


2hk% 


The non-oscillating part of (14) is given by 


ff( 0 ) = 

^ ZZ 


UeC^VF 

hkp 


1 - 


5 ^F 


3C| 


(14) 


(15) 


(16) 


1 /t ioj + Tj ^ 

where rie = kp/{3TT'^) is the total density of electrons per 
chirality. And, for the oscillating part we have 

3/2 

, (17) 


(/) ^ nee VF 

^ ZZ 


1 


XI 


hkF iw-b 1 /t 27r sinh A/ \ I J 

where A = tt'^T/{ hkFV fCf)■ In the DC limit and in the 
absence of magnetic field, C,f = 0, Eqs. are re¬ 

duced to a standard Drude formula appropriately modi¬ 
fied for Dirac spectrum: 


(To = 


UeC^T 

hkF/vF 


(18) 


In finite magnetic field the second term of (16) de¬ 


scribes an ideal conductivity. In the absence of chirality 
flipping this conductivity diverges in static limit w —> 0 . 


In more realistic models the process of chirality flipping 
are always present and one should replace rj —)■ l/r«, 
where Ty is a mean chirality lifetime. As the scattering 
with and without changes of chirality are due to very dif¬ 
ferent processes one should expect the ratio r^/r to be 
significant. Both r and Ty can in principle, be extracted 
from optical conductivity measurements. 

There are two small parameters in the regime of in¬ 
terest of this work. One is (f, be., the weakness of the 
magnetic field compared to the Eermi scale. The other is 
the smallness of temperature compared to the Eermi en¬ 
ergy. We do not, however, make any assumptions on the 


relative size A of these small parameters. In deriving (14 


17) we kept the leading (H-independent) and next to the 


leading terms of the expansion in (^f but restricted the 
expansion only to the leading term in TjsF- This is why 
the only temperature dependence in ( MpT ) is through 
the parameter A. This means that we omitted all correc¬ 
tions proportional to T/ef which could be comparable 
to the ones proportional to (^f- The former corrections, 
however, are not universal and do not affect the magnetic 
field dependence of the conductivity. 

A very convenient way to exclude the non-universal 
temperature corrections is to study the ratio a^z^B)/ uq. 
In DC limit (w —>■ 0), it is given by 

(^zz{B) , o f 4 \ 2 3 


o-Q 


= l-b3 


T I Cj’ 4 

TO/ TT 


(19) 




s—^dI 


XI 


sinh XI 


" t ) 


^/2 


COS 


ttI 

2 C,f 


In the last equation we introduced the Dingle factor 
Xf) = ttT/ ihvFkFCF), which accounts for the smearing of 
LLs due to impurities. In the case homogeneous sample 
r = K/tq with the “quantum time” tq determined by 
impurity scattering and equal to quasiparticles’ lifetime. 

If either A 1 or Ad » 1, i.e., the temperature or 
smearing of Landau levels is larger than the gap between 


Landau levels, the oscillations in (19) disappear and the 


conductivity is given by the first line in (19). Eor smaller 


temperatures and Landau level smearing oscillations ap¬ 
pear and become less and less harmonic with a further 
decrease of both A and Ad- 

In Fig. we plot the magnetoresisitivity given by the 


inverse of expression in Eq. (19) for parameters consis¬ 


tent with the recent experiment on Cd 3 As 2 m- Com¬ 
parison with the Fig. 4B of m shows that the approach 
to magnetotransport in Dirac semimetals developed here 
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FIG. 2: Longitudinal magnetoresistance as a function of the 
magnetic field. We used numerical values consistent with m- 
kp = 3.8 X vf = 9.3 x 10®m/s, r = 8 x T = 

2.5K. The plots are made for three values of tq/t — 0,1,16 
and for r„ = lOr. 

describes qualitatively the emergence of quantum SdH 
oscillations and the tendency to negative magnetoresis¬ 
tance at strong magnetic fields (but still small (p) ob¬ 
served experimentally in Cd 3 As 2 However the 

Cd 3 As 2 data exhibit also a strong positive magnetore¬ 
sistance present in weak magnetic fields m The more 
thorough comparison of our theory with experimental 
data requires an explanation of the positive magnetore¬ 
sistance. Since a (much weaker) positive MR has also 
been observed in the weak magnetic field region in Dirac 
semimetals ZrTes and NasBi, it is likely that this effect 
is generic for three-dimensional chiral materials. Possi¬ 
ble explanations include the magnetic field dependence of 
impurity screening, the weak antilocalization and surface 
effects resulting from Fermi arcs m- We leave the sys¬ 
tematic treatment of these effects, as well as the study of 
microscopic mechanisms of chirality relaxation, for future 
studies. 
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The negative MR was not observed in m- 
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APPENDIX A: BOLTZMANN EQUATION 


In this section, we will derive the expression (11) by 


solving the stationary Boltzman equation. Rewriting Eq. 
in spherical coordinates and plugging the formula for 


scattering rate ( 10 ) into it, we obtain: 


d 

iuiT + 2kC,kVk-^ ) g{k) - Vkk = 
3T 


(A.l) 


= VkXCkZ + z^ I d^k'T'ig' - g ) ^ - s')- 


(1 + fc' • k) 


IGtt^ 


v{e)T{e) 


Integrating Eq. (A.l) over the solid angle, we obtain: 

d-nxCkZ 


/ d0d(cos0) T(A;)g(fc) = 

/g2 iw 


(A.2) 


It is obvious from here that there are no stationary so¬ 
lutions for Eq. 0 when w —>■ 0. The azimuthal symme¬ 


try allows to find solutions of (A.l) that are independent 


of (j) for gz- After the integration over (k', (/>'), we end up 
with: 

iujTg^ = Vkicos6 + xCk) + J d{cos0')T'{g'^ - g^)T 
X ^( 1 -I-COS 0 COS 0 '). (A. 3 ) 

The easiest way to solve this equation is to expand Tgz 
in terms of Legendre polynomials and use their orthogo¬ 
nality conditions. Thus, 


Tg, = ^( 2 Z + l)ai(fc)Pz( COS 0 ), 
1=0 

XCkVk 


lU) 


cl , ( 1 -Cfc^) 

ioj iuj-\-T~^ 


VkCos9, (A.4) 


where ag is obtained through (A.2|. 


APPENDIX B: DISCRETENESS OF LANDAU 
LEVELS 

Quantum effects in the conductivity can be imple¬ 
mented through the Bohr-Sommerfeld quantization con¬ 
dition. The prescription here is the same one used in the 
old quantum theory; given a classical system, we intro¬ 
duce quantum effects by imposing that canonical vari¬ 
ables satisfy: 

j) Pi dqi = 2TTh {tii + jindy), 


where 7 is a curve in phase space in which the Hamilto¬ 
nian is a constant and indy is the Maslov index of y. 

However, in the presence of a nonvanishing Berry cur¬ 
vature, the perpendicular components of k^^ fail to be 
canonically conjugated. Instead, 

1 / p \ p 

+ = (B.l) 

Following the same recipe and using that B = Bz, the 
discreteness of Landau levels can be imposed by assuming 
that: 


i( 


1 -H Q JB 
h 


2tt 


k X dk = ^{n+ jindy). (B.2) 


It implies the area quantization - in units of 27r/£^ 
- for the section of the Brillouin zone with kz constant. 
We can find the surfaces with constant n in fc-space by 


solving equation (B.2): 


1 . 

4^ 


If we impose that n = 0 is the smallest possible integer 



. to 

1 

(B.3) 

2 

V + ^1/' 

I 

4Cfc 

(l — cos^ 9 — 2xCk cos 9 ) . 

(B.4) 


solution of (B.4) and use the fact that indy € Z; the only 


possible values of the Maslov index are {—2, —1, 0}. 


APPENDIX C: SDH OSCILLATIONS 


In this section, we will apply the Bohr-Sommerfeld 
quantization prescription to introduce quantum effects in 
the conductivity. Assuming that the only contribution to 
transport comes from the discrete levels, the conductivity 
per chirality becomes: 


9 00 

cjA) = -^ 

27r2 

n= —00 


^ y d(cos6l)A:2^(cos6i-f xCfe) 

= — 00 _ 

X S[n- I (sin^6>/Cfc - 2xcos6> - md'y)]vkgz, 


(C.l) 


where gz is given in (II). We have used that there is 
no real solution for (B.4) when n G Z_. Therefore, all 
surfaces for negative integer n are outside of the integra¬ 
tion range. The integral over k is performed near Fermi 
surface. 


M With respect to B. 



























Using the Poisson formula, 


where 


6 {x-n)= Y 


A27Tln 


1 — — C 


the conductivity can be rewritten as: 


O OO 

^(x) — 

27r2 

i=- 


E 


X 


(0 


r{l) 


+ - 


iuj + Tv ^ ' ioj + T 1 


(C.2) 


(C.3) 


where, 

Xf ^ =Jdk k^l ^ xCk 

1 

X J d(cos6*)(cos0 + xCfc) > 

-1 

Xf =J dkk'^vl ^ (1 - Cfe)e*T(U-^+C'=-”d7) (C.4) 

X /d(cos0)^2^^i^ cos0e-*TUo>^e+xCU= 


-1 


1 + xCfeCos 6 > 


The integral in Eq. (C.3) accounts for the intervalley 


scattering and can be easily calculated: 


i+xCfc 

Qm = J 

-1+xCfc 

Solving for odd values of m: 


1 m I ■ i 

axx exp [ —%-—x 
2 Cfc 


(C.8) 


Qi = ^ sin(7r;x) = 0 , 

TTL 


23 = 


4x*(-l)'Cfc 


ttI 


• e 2 


(Cfc +Cfc) 


However, Cfc23 = C’(Cfc) and such term can be ne¬ 
glected. Let us now focus on m even. They can all be 
obtained through Qo as follows: 


Qm = 




\ TT J 


dlm/2 


Qo ) 


where I is set to be a non-zero integer at the end of the 
calculation. Clearly, Q 4 = O(C^) and we only need to 
calculate Q 2 . The integral Qo can be written as: 

i-i-xCfc i-xCfc 


Qo= / dxe~'^"Y 


dcce . (C.9) 


Xf ^ j ^ 1,0 ■ 

In the equation above, we have used that 

S(k — kp) 


(C.5) 


dh 

ds 


hvp 


The terms with I ^ 0 vanish since the chirality relax¬ 
ation mechanism that we have considered only accounts 
for the scattering between the zero-modes. 

Let us now consider the contribution for the intravalley 


scattering coming from Eq. (C.4). For I = 0 the inte¬ 


gral can be performed analytically, however, we are only 
interested in the range where the semiclassical picture is 
valid. If we restrict ourselves terms up to 0{Cp), we end 
up with: 


4°' = 




In order to calculate T 2 '’ for / ^ 0, it is convenient to 
define x = cos 6 -|- xCk ■ Expanding the integrand up to 
C>(Cfc), we find that: 


(C.6) 


-rlO /" u; ; 2 2 ^/o 
1^2 = i dkk v^— exp 


7tI 

iCk - indy - I/Cfe) 


X [—xCfcQi + (1 + Cfe)Q2 — CfcxQs + CkQi]^ (QX) 


Let US focus on the right hand side of Eq. (C.9). Thus, 
l±xa 

■). 


1 . /^Cfe 1 , /2Cfe i:, f ii-|(l±xCfcP 


= 2 VT- 2 Vi?"'»(" 


2Cfe 


where we have defined: 


,-”i-l/2g-iy 


/ OO 

dyy~ 

After integration by parts, one can show that: 

Prn+l{t). 

Since we are restricting ourselves to terms up to 0{C,1), 
we obtain: 


Qo = 


Kk , 2 a(-l)'e-*T(U-^+a) 



-h 


tI 






Q 2 =—v^ + 

TT \ 


7tI 


e-*¥(U7^+a) {i+‘^ 
TTl 


Plugging all determined values for into (C.7): 
(0 _ _ ‘^kpVpCF ^i-Kljl-kind-y) ( ■ , 


-^2 ~ 


HttI 

y2e*f(5-“d7/) 

h'Kp/'^ 


dek VkC, 


2 ,„.X/ 2 y 2 ei¥(c.-‘+c.) 


de 
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Here we have used that we can invert the dispersion 
relation and write k{s). The energy integral is performed 
at the vicinity of the Fermi surface. Since we assume 
that T/ep ^ 1, all the integrand besides the oscillating 
exponential is consider to vary slowly in the temperature 
range. In addition to that, we must expand the exponent 
near the Fermi energy. Keeping only linear deviations in 
the exponent, we are left with: 


t (0 

-^2 


HttI \ irl J 

if [i+i(Cp^-ind7)] / \ T 

-2fa-(-rj J "WTF’ 


In the Eq. ( C.1I| ), we have defined t = (e — /r)/T and 
A = tt'^T/{ hvFkpCF)- The integral can be solved using 
the residue theorem, and its value is given by: 


OO 



— OO 


g(l+iAi/'n-)t 

(e‘ + 1)2 


dt = 


a; 

sinh \l 


Therefore, the conductivity can be expressed as: 




(0 


cos 


1=1 


ttI f \ 
Cf 


— indy ) + 


4_ 
(C.12) 


where 


( 0 ) _ n.ee vp 
hkp 


3C| 

iuj + r. 


-1 


1 

iu) + T~ 


(C.13) 


X ( 1 + 2^-^ ^i+ind7,o —^Cf + ^(ind7)^CF 


and 


_ Uee^vp 1 5 XI / 2Cf A 

kkp iu! + l/r 2 p sinh Al \ I J 


3/2 


. (C.14) 


In the Eq. (C.13), we have used that: 




Although indy can in principle be obtained by the 
WKB calculation, we assume indy = 0 in the main text. 


APPENDIX D: DINGLE FACTOR 


that the density of states have sharp peaks at each Lan¬ 
dau level. However, this is not true in a more realistic 
scenario. The presence of impurities breaks the energy 
degeneracy of the Landau levels and as a net result they 
get smeared by the presence of impurities. 

The assumption that kp£g/{vpT) ^ 1, allows us to 
disregard corrections to the plane-wave scattering due to 
the magnetic field^^. Within this approximation, the den¬ 
sity of states is still isotropic, however, it gets a contribu¬ 
tion coming from the smearing of energy levels, namely: 


f d 3 fc 3 [G«(e,fe)] , 

J ]BZ 

7(0 = i/ A. 


(D.l) 


[^_e(fc)] 2 +r 2 (^,A:)- 

In the limit when F —> 0, we recover the well-know 
result 

47rfc2 


KC) = / <5(e - e{k)) = 

J MZi 


hvk 


e=5 


In fact, the smearing of the energy levels can be intro¬ 
duced by following replacement: 

- eik)) ^ 1 

TT - e{k)Y +T^{^,k) 

We can thus rewrite Eq. ( C.10| ) in a more convenient 
way: 


4^ = -^k,^i-ly[^ + 

nTTL 


2Cp\ ^f2e^l 


ttI 


7r^3/2 


(D.2) 


X / d^ / dkS{^-eik))evlct^^^ 






Here, we have set indy to zero in order to shorten up 
the notation since this factor brings no extra difficulties. 
The choice of the lower limit of integration is for later 
convenience. As previously mentioned, the smearing can 
be taken into account by replacing the delta function in 
the integral above by a Lorentzian distribution. There¬ 
fore, let us focus on: 


SiiO 



dk 


r(^,A:)e*"2 (Cfc'+a) 3/2 

K-^(fc)]2 + r2(e,fc) 


(D.3) 


In the treatment of quantum oscillations, the Dingle 
factor in Eq. (20) comes from the smearing of LLs due to 
impurity scattering. In the previous section, we assumed 


Otherwise, we must consider the whole matrix elements of the 
impurity potential in the presence of magnetic field. 
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The cutoff 1/(\/2 £_b) guarantees the integral conver¬ 


gence. One can solve Eq. (D.3) using the steepest descent 


approximation. For that, let us analytically continue the 
integrand and define = Cfe Hence, 


SiiO = 




dzH 


V^Ib) 




where 


m,k) = r, 


m,k)vi 


[i-s{k)]^+T^i,k)- 

Expanding the exponent near z = 1, we obtain: 


(-lY 




(D.4) 







The contour C is defined by 3fi[(z — 1)^] = 0 together 
with 3[(2: — 1)^] > 0 and \z\ > 1. Let us assume for 
simplicity that r{^,k) = r(^). Using that only C ~ M 
contributes to we find that: 




(-l)'uo r kYvkCk^ _ ( i'Kl\ 






where 


and k is taken to be k{^ + iT). Plugging it into 
end up with: 


(0 _ I l'^ ^ ‘2 Cf\ I 


IX’ = - 
^ HttI 




7v‘^Pep 




2h7T 


-kpVp 


"t) 


3/2 


dt 


(e* -I- 1)^ 


Here, we have defined Xd = ttT/{ hkpVFCF) and ne¬ 
glected terms of However, from Q, 

T/sf T/{hvFkF) ^ Cf 

and consequently (^y^T/e <C Cf- Therefore, the second 
term in 1^ can also be neglected within our approxima¬ 
tion. 

The only modification in the conductivity expression 


coming from the smearing of LLs occurs in Eq. (D.6|, 
which must be replaced by: 


gd) - 3 XI /2CFy^^^ 

hkF iw-|-1 /t 2f sinh AZ \ I J 


Vo = Vk\, _ 

V2£b 






















